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Q I Abstract 

We propose a sufficient condition for invertibility of a polynomial 
mapping function defined on a cube or simplex. This condition is ap- 
plicable to finite element analysis using curved meshes. The sufficient 
' condition is based on an analysis of the Bernstein-Bezier form of the 

QQ . columns of the derivative. 

o 

cn 

5 ■ 1 Invertibility of polynomial mapping func- 

P. ; tions 

> . 

^ . In finite element analysis, it is common to subdivide the domain into elements 

^ ! that are images of a reference domain under polynomial functions. This 

approach gives rise to the popular isoparametric elements [6]. The reference 
domain is the unit square P = : < ^,?7 < 1} or triangle = 

{i^,v) : < ^,77;^ + r/ < 1} in R2 or the unit cube P = {i^,r]X) ■ < 
^,V,C< 1} or tetrahedron = {(e, r/, C) : < e, r?, C; ^ + ^ + C < 1} m R'- 
Polynomials defined on A'', d = 2,3, generally include monomial terms 
up to degree p for some p > 0. On the other hand, for the unit cube I'^, the 
monomial terms are generally up to degree p individually in each coordinate. 
Therefore, for the rest of the paper we use p to denote the maximum total 
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degree of polynomials on or A^, and the maximum degree in individual 
coordinates of polynomials defined on P or I^. 

Let F : U R"' be a polynomial function, where U is either J"' or 
A"^ and o? is 2 or 3. Suppose F is injective on [/, implying that it has an 
inverse G defined on F{U) such that G(F(u)) — u. If such a G exists 
and is smooth, then we say that F is invertible. Invertibility is important 
for correctness of a finite element mesh. This definition of invertibility is 
sometimes called "global invertibility" in order to distinguish it from "local 
invertibility." The function F is locally invertible if its derivative, denoted 
as J, is nonsingular on the entire reference element U. Local invertibility 
is necessary for global invertibility but is not sufficient. A function can be 
locally but not globally invertible if F{U) "wraps around" (e.g., consider 
F{r,9) = r{cos9,sm9) for (r, 6*) G [1,2] x [0,47r]). Additional sufficient 
conditions for global invertibility in the general setting are considered by 
Ivanenko [5]. Beyond local invertibility, one would also like an upper bound 
on the "condition number" 

cond(F) = max||J(u)|| -maxllJlu)"! (1) 

of F. Even more strongly, one would like bounds on the higher derivatives 
of F in terms of cond(F) to satisfy the Ciarlct-Raviart [2] conditions on 

convergence of order-p finite element approximation. 

Unfortunately, even the first problem, namely, determining nonsingularity 
of the Jacobian, is a difficult problem. We know of simple necessary and 
sufficient conditions only for the following special cases: 

1. linear polynomials, that is, p = 1, defined on A*^ of all dimensions, 

2. quadratic polynomials on A^ in the case = 2, in the special case that 
F is linear on two of the three edges of the reference triangle [6] , and 

3. bilinear elements (i.e., p — 1) defined on P [14]. 

For all other cases, we know of only separate necessary and sufficient 
conditions, and this paper also proposes only a sufficient condition. To our 
knowledge, the only necessary condition for the general case proposed in 
the literature is that det(J) have the same sign (strictly positive or strictly 
negative) at some finite list of test points. In engineering applications, it is 
common to test invertibility at the Gauss points used for quadrature on the 
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element [13]. This test is known not to be sufficient; see further remarks on 
this point below. 

Sufficient conditions have been proposed for a few settings. For example 
Ushakova [15] presents a sufficient condition for the case p = 1 on domain 
I^. Her sufficient condition is interpreted as requiring that certain tetrahedra 
formed by choosing subsets of four of the eight vertices of F{P) have positive 
volume. Sufficient conditions for quadratic triangles and tetrahedra have 
been proposed by Salem, Canann and Saigal [9, 11, 10, 12]. 

In this work, we present a sufficient condition to ensure both local and 
global invertibihty of F for any of the four domains I'^, A'^, d = 2,3 and for 
any polynomial degree p. Our condition appears to be weaker (i.e., not able 
to certify invertibihty in more cases) than existing sufficient tests for specific 
reference domains and values of p but is considerably more general. 

Our condition is based on writing J in Bernstein-Bezier form and then 
considering the convex hull of the derivatives at the control points. In the 
next section we describe the new sufficient condition and establish that it is 
sufficient for local invertibihty. In Section 3 we provide an equivalent charac- 
terization of the sufficient condition that is amenable to efficient testing. In 
Section 4 we use the equivalent characterization to prove that the condition 
also implies global invertibihty. Our condition has a desirable property that 
we term "affinity," which we define in the last section. 

2 Bernstein-Bezier form 

Bernstein-Bezier (BB) form is a popular way to write polynomials in computer- 
aided geometric design [4]. A univariate polynomial of degree p in BB form 
would be written: 

nO = E/i(l-erT.yT^ (2) 

and has natural parametric domain ^ e [0,1]- A bivariate polynomial with 
maximum degree p individually in ^, 77 is written in the form: 
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and has as its natural parametric domain 77) e 7^. A trivariate polynomial 
with maximum degree p individually is written 

(4) 

and has as its natural parametric domain {^,ri,() G 

A bivariate polynomial with total degree at most p is written in the BB 
form 

n«..)^|:|«V(i-^-,r-^^^;;^ (5) 

and has natural parametric domain (C^v) ^ 

Finally, a trivariate polynomial with total degree at most p is written in 
the BB form 

^' = Sirio '-'^'"'^''^ 

(6) 

and has natural parametric domain {^.r/.Q G A^. 

In all five cases, the vectors (in ID), fjj G (in 2D) or fjj^fc G (in 
3D) are called the control points. A fundamental theorem about BB form is: 

Theorem 1 Let U be the natural parametric domain of BB form for the five 
cases listed above. Then F{U) is contained in the convex hull of the control 
points. 

Farin [4] proves this as a consequence of the deCasteljau algorithm for 
evaluating F, but here is a sketch of a more direct proof. One observes that 
on the natural parametric domain, all the factors in the summations like ^, 
(1 — ^ — r/), etc., are nonnegative. Furthermore, one observes that the value 
of all the sums, if the control points are excluded, is exactly 1. For example, 
consider (2) without control points: 

E(i-erTTT/-T7 

i\{p-i)\ 

This sum is identically 1, as seen by considering a binomial expansion of 
(^ + (1 — ^)Y. Similar argument apply to (3)~(6). Thus, in the natural 
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parametric domains, the above summations may be regarded as weighted 
averages of the control points, proving the theorem. 

The next feature of BB form is that once a function is in BB form, the 
derivatives can be easily put into BB form. In all cases, the derivative is a 
BB expansion of one degree lower, multiplied by p, and with control points 
that are finite differences of the control points for F in the direction of the 
variable being differentiated. Thus, for example, in the case of (5), we have 

QP p-lp-j-l (v-lV 

and 

In this case, the list of p{p + l)/2 vectors of the form p(fj_|_i j — fjj) are control 
points of dF/ d^, and analogously ]9(fjj+i — fjj) are control points for dF/drj. 
Let us denote these two lists of control points and respectively. Similar 
expressions hold for the other four BB forms described above. In 3D there 
is a third list G^. We now state our condition: 

Condition 1 In the case d = 2, the matrix [u, v] is invertible for every 
u e hull(G^) and every v e hul^G^). In the case d — 3, the matrix [u, v, w] 
is invertible for every u e hul^G^), v e hul^G^,), w e hul^G^). 

In this condition, "hull" denotes the convex hull. The main result of this 
section is as follows. 

Theorem 2 // Condition 1 holds, then the matrix J is invertible on the 
entire reference element. 

The proof of this theorem follows from the arguments in the previous 
paragraphs: Condition 1 is sufficient since the actual Jacobians that occur 
on the domain have their first columns chosen from hull(G^), etc. 
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3 A computational characterization of the suf- 
ficient condition 



Condition 1 specifies our sufficient condition in somewhat nonconstructive 
terms. In this section we provide a computational means to verify the suffi- 
cient condition. 

Theorem 3 In the case of R^, Condition 1 is equivalent to the following 
condition: 

• there exists a vector h e such that for all f e U G^, h-^f > 0, 
and 

• there exists a vector h e R^ such that for all f & G^, h^f > and for 
all f e Gr,, h^f < 0. 

Proof. Let u, v be arbitrary in hul^G^) and hull(Gr,) respectively. This is 

equivalent to saying there exist nonnegative summing to 1 such 

that u = Q;ifi + - • ■ + Q;„f„, where fi, . . . , f„ is an enumeration of G^. Similarly, 
there exist nonnegative (3i, . . . , {3^ summing to 1 such that v = /5igi + • ■ ■ + 
PmSm, where gi, . . . , is an enumeration of G^. The condition that [u, v] is 
invertible is equivalent to saying that u, v are independent, i.e., that there do 
not exist S, 7, at least one nonzero such that Su + 7V = 0. Suppose they are 
dependent, and let 6, 7 be the two coefficients of dependence. Without loss of 
generality, 6 > 0. There are now two cases: either 7>0or7<0. If7>0, 
define ai, • • • , to be Sai, . . . , 5an and define - ■ ■ ,/3m to be 7/3i, • • • , 7/3„. 
Then the condition that the a's are nonnegative and the assumption that 
5 > is equivalent to the hypothesis that ai, . . are nonnegative (with 
no restriction on their sum). Similarly, the condition on the /S^'s is that they 
are nonnegative. 

Thus, in the case that 7 > 0, a dependence is equivalent to the existence 
of nonnegative ctj's and PiS, not all zeros, such that 

n m 
i=l i=l 

This problem is solved by linear programming. By Farkas' lemma [16], there 
is a nonnegative solution to this problem, with not all the coefficients zero, if 
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and only if there does not exist a vector h such that h^f, > and h^gj > 
for all fj's and gj's. 

Now wc turn to the case of dependence when 7 < 0. In this case, using 
the analogous argument, there is no dependence provided that there is a 
vector h such that h^fj > for all i and h^(— gj) > for all i. I 

This shows that in the case of R^, Condition 1 can be tested by solving 
two linear programming problems over to find h and h. In fact, there 
is a much simpler algorithm, namely, compute the args of all the points in 

U Gr,. To determine whether h exists, one checks whether the max arg 
differs by less than tt from the min arg. (The branch cut for defining the 
arg function must be chosen outside the min- max range). A similar test can 
determine whether h exists. 

In the three-dimensional setting, the algorithm is not so simple but it is 
still linear time. Following the same approach as in the previous proof, we 
see that in 3D, Condition 1 is equivalent to the following conditions: 

• there exists a vector hi e R^ such that hf f > for any f e G^UG^UG^, 

• there exists a vector h2 G R^ such that li^f > for any f e G^ U G^ U 
(— G^), where — G^ means {— f : f G G^}, 

• there exists a vector h.3 G R^ such that h^^f > for any f G G^ U 

{-Gr,) U G^, and 

• there exists a vector h4 G R^ such that hjf > for any f G G^ U 
(-G,)U(-Gc). 

Each of these can be verified with linear programming, which is linear time 
in three dimensions [8]. Alternatively, they can be verified in 0{n\ogn) time 
using a convex hull algorithm [3] . 

4 Global invert ibility 

We are now in a position to prove that Condition 1 in fact implies global 
invertibility of the mapping function. To demonstrate global invertibility, 
we will show that F is injective, i.e., for all u, v in the reference domain, 
u 7^ V ^ F{u) 7^ F{v). Injectivity is sufficient for global invertibility: if 
an injective smooth function has a nonsingular derivative at all points of 
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its domain, then the global inverse is also smooth by the inverse mapping 
theorem. The nonsingularity of the derivative is already established. 

Theorem 4 If Condition 1 holds, then F is globally invertible on the refer- 
ence element. 

Proof. Let us start with the proof of injectivity in two dimensions. Choose 
u, V in the reference element such that u 7^ v. Let w = v — u. Then by the 
definition of directional derivative, 



Assume without loss of generality that Wi > 0, where {wi,W2) denotes the 
entries of w. (The case Wi < is handled by exchanging the roles of u, v.) 
Now there are two cases: either ■w;2 > or ■w;2 < 0. Suppose first that W2 > 0. 
Because we assume Condition 1 and because we have established Theorem 3, 
we conclude that there exists a vector h e such that h^f > for all 
f e U Gri- By convexity, this implies h'^f > for all f G huII(G'^ U G,,). 
This implies that h^J{^,r]) is a vector both of whose entries are positive 
for any ^, in the reference domain since both columns of J are taken from 
hull(G^ U Gjj). Since w has both entries nonnegative, and at least one entry 
of w is positive (because u 7^ v), this means li^ J{^,ri)w > for all $,,1] in 
the reference domain. Thus, taking the inner product of both sides of (7) 
with h shows that h^(F(v) — F{u)) > 0, and in particular, F{v) ^ -F(u). 

The second case is W2 < 0. The argument is analogous, except we use h 
instead of h. 

Finally, the extension to 3D follows the same lines. We assume w\ > 0. 
Depending on the signs of entries W2, W3 of w, we select one of hi, h2, ha, h4. 
I 

5 Discussion 

We say that a sufficient condition for invertibihty has the "affinity" property 
provided that it is always satisfied if the mapping function is a nondegenerate 
affine linear mapping, or a sufficiently small perturbation of a nondegenerate 
affine linear mapping. 
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Affinity is a desirable property since the affine linear mapping is obviously 
the easiest case for invertibility. Furthermore, if a sufficient condition for in- 
vertibility has the affinity property, then it can be turned into a necessary 
and sufficient condition if applied to each subcell of a sufficiently fine sub- 
division of the original domain. This is because an invertible mapping on a 
fine subdivision will behave like an affine mapping plus a small perturbation 
on each subcell, and the small perturbation tends to zero as the subdivision 
gets finer. 

It is fairly easy to see that our test has the affinity property. Recall 
that the BB control points for a constant function are all identical. Thus, 
if F is affine linear, then G^, and are all singleton sets. Assuming 
that F is affine and invertible means that the condition must be satisfied. 
Furthermore, a sufficiently small perturbation of F will give G^, G^, and Gq 
a positive radius, but the condition of separability by planes will still hold. 

Ushakova's condition also has the affinity property. We suggest that any 
proposed sufficient condition should possess the affinity property in order to 
be considered useful. 

Ushakova's sufficient condition for the p = 1 and U = case is less 
restrictive (i.e., is able to validate more elements) than the condition proposed 
here. 

We conclude with a few open questions raised by this work. 

1. In the case p = 1 and U = I^, our condition is more restrictive than 
Ushakova's. On the other hand, it is not obvious how to generalize 
her condition to higher degrees or to tetrahedral elements. It would 
be interesting to combine her techniques with the ones in this paper to 
come up with less restrictive conditions for higher degrees. Ushakova 
(private communication) found that experiments with randomly gen- 
erated elements indicate that for p = 1, U — the conditions stated 
here overly restrictive in the sense that most of the valid (invertible) 
elements among a randomly generated test set would not satisfy the 
sufficient condition proposed here. 

2. In practical application of isoparametric mesh generation, one defines 
the polynomial function only on one face of the boundary of the ele- 
ment (namely, the face adjacent to a curved exterior boundary). The 
remaining coefficients defining F can be determined by a formula such 
as the formula proposed by Lenoir [7], which has certain theoretical 
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guarantees. How does the sufficient condition here speciahze if we as- 
sume that interior degrees of freedom of F are determined by Lenoir's 
formula? Our very preliminary computational tests seem to indicate 
that using Lenoir's formula seems to make the element more amenable 
to our sufficient condition (i.e., it seems to minimize the diameters of 
the sets G^,G^, compared to other choices for interior degrees of 
freedom) . 

3. The test proposed here could be strengthened to obtain an upper bound 
on cond(F) defined by (1) in terms of the numerical values of h-^fj as 
fj ranges over the sets G^, G^, G^ and h is the vector defining one of 
the halfspaces in Theorem 3. This raises the possibility of a sufficient 
condition not only for invertibility but also for confirming the Ciarlet- 
Raviart conditions. Recall that the Ciarlet-Raviart conditions require 
that higher derivatives of F are bounded in terms of cond(F) . 

Note that other condition numbers besides cond(F) may be useful in 
practice. For example, Branets and Garanzha [1] have developed a 
distortion measure related to structured grid generation. It would be 
interesting to get bounds on all of these condition numbers in terms of 
the the numerical values of h^fj. 
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